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ABSTRACT 

We report the discovery of circularly polarised radio emission from the radio-jet X-ray binary SS 433 
with the Australia Telescope Compact Array. The flux density spectrum of the circular polarization, 
clearly detected at four frequencies between 1-9 GHz, is of the form V oc jy-°- 9 ± 01 . Multiple components 
in the source and a lack of very high spatial resolution do not allow a unique determination of the origin 
of the circular polarization, nor of the spectrum of fractional polarization. However, we argue that the 
emission is likely to arise in the inner regions of the binary, possibly via propagation-induced conversion 
of linear to circular polarization, and the fractional circular polarization of these regions may be as high 
as 10%. Observations such as these have the potential to investigate the composition, whether pairs or 
baryonic, of the ejecta from X-ray binaries. 

Subject headings: Radio continuum:stars - Stars:individual:SS 433 
accretion, accretion discs - ISM: jets and outflows 



1. INTRODUCTION 

High-velocity synchrotron-emitting jets are commonly 
observed from Active Galactic Nuclei (AGN; e.g. Os- 
trowski et al. 1997; Shields 1999), and X-ray binary sys- 
tems (XRBs) containing both black holes and neutron 
stars (e.g. Hjellming & Han 1995; Fender 1999 and refer- 
ences therein). The composition of the jet plasma, whether 
electron-proton (e~p + ) or electron-positron (e~e + ) re- 
mains a fundamental yet unanswered question in nearly 
all cases. 

SS 433 is one of the most celebrated of Galactic ob- 
jects. The source is an XRB consisting, most probably, 
of a mass-losing star in a 13-day orbit with a stellar- 
mass black hole or neutron star. The system produces 
bright quasi-continuous radio jets which precess with a 
period of ~ 162.5 days (Vermeulen 1989; Vermeulen 1992; 
Brinkmann 1998). Moving optical emission lines (Mar- 
gon 1984 and references therein) indicate a jet velocity of 
(3 = v/c = 0.26, confirmed by both VLA and VLBI ra- 
dio observations (Vermeulen 1992 and references therein). 
These optical lines, and their X-ray counterparts (Kotani 
et al. 1996 and references therein), are the only direct ev- 



idence for the existence of baryonic material (i.e. e~p + ) 
in a jet from any X-ray binary. 

Recent progress towards determining the composition of 
the plasma in jets from AGN has been made by the detec- 
tion and modelling of a circularly polarised radio compo- 
nent from the quasar 3C 279 (Wardle et al. 1998). Wilson 
& Weiler (1997) argue that radio circular polarization up- 
per limits for the Crab supernova remnant come close to 
determining the positron content of the nebula. In ad- 
dition Bower, Falcke & Backer (1999) and Sault & Mac- 
quart (1999) have recently detected circularly polarised 
radio emission from Sgr A* at the Galactic Centre. 

In this paper we report the detection of circularly po- 



larised radio emission from SS 433, the first from any XRB, 
at four radio frequencies. This observation has the poten- 
tial to be the benchmark against which other jets from 
other XRBs may be compared in an effort to determine 
whether they produce e _ p + or e _ e + jets. 



2. OBSERVATIONS AND RESULTS 



The Australia Telescope Compact Array (ATCA) con- 
sists of six 22 m alt-az antennas near Narrabri, New South 
Wales (Frater, Brooks, & Whiteoak 1992). Each ATCA 
antenna is equipped with two wide-band feed horns, and 
each feed horn is equipped with two pairs of orthogonal- 
linear probes. This allows both orthogonal polarizations 
at two separate frequencies to be observed with each feed 
horn at the same time. Observations were made using 
the ATCA 'continuum-mode', which gives bandwidths of 
128 MHz simultaneously at each of two frequencies, and 
four correlation products (XX, YY, XY, YX). The obser- 
vations on 1999 May 10 were made centered on 4.80 GHz 
and 8.64 GHz. On 1999 May 20, observations at 1.38 GHz 
and 2.50 GHz were alternated with those at 4.80 GHz and 
8.64 GHz, with a cycle time of ~ 25 minutes. Calibra- 
tion sources, selected from the ATCA Calibrator Source 
Catalogue (Reynolds 1997), were observed every cycle. 
The ATCA primary calibrator PKS 1934-638 was used to 
calibrate the bandpass and to set the absolute flux scale 
(Reynolds 1994). 

Data reduction was performed with the Miriad package 
(Sault, Teuben, & Wright 1995). Prior to gain calibration, 
the xy-phase correction measured by a noise-diode sys- 
tem was applied, and the data were corrected for a small 
(otherwise unmodelled) field rotation due to the antenna's 
pointing model (Kesteven 1997). The main calibration 
step involves simultaneously solving for time-dependent 
complex-gains, time-independent residual xy-phase, and 
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Fig. 1. — Fractional circular polarization for SS433 and calibrators for 1999 May 20 observations. The scale is set by assuming PKS 1934-638 
has Stokes V=0 at all frequencies (see text). 




Fig. 2. — Observations of SS 433 in 1999 May. Top panel illustrates GBI (quasi-)daily monitoring at 2.3 & 8.3 GHz; our best estimates 
of the ejection dates corresponding to the two major flares are indicated by arrows. The lower panels indicate mean flux densities in total 
intensity, linear and circular polarizations for the two epochs of ATCA observations. 
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Fig. 3. — Progression of internal and external absorption and depolarization effects as ejecta move away from SS 433 (Hjellming & Johnston 
1981; Vermeulen 1992; Paragi et al. 1999). The circles at the top and bottom of the schematic indicate our best estimates of the physical 
locations of the components responsible for the radio flares which dominate the total flux density, on 1999 May 10 and May 20 respectively. 



time-independent polarization leakages for each feed, as 
well as the linear and circular polarization of the calibra- 
tion source. A good parallactic-angle coverage is required 
in order that the leakages and calibrator linear polariza- 
tion can be decoupled in the solution process through the 
relative rotation of the feeds and parallactic angle (Con- 
way & Kronberg 1969). 

Circular polarization data requires calibration using the 
'strongly-polarized' equations (Sault, Killeen, & Kesteven 
1991), which include terms in leakage x (Stokes Q,U). 
To obtain such a solution requires a calibration source 
which has a few percent linear polarization (for this ex- 
periment PKS 1908-202) so that there is sufficient signal 
in the second-order terms. The accuracy of the leakage 
solution was estimated by repeating the entire calibration 
procedure for the 1999 May 20 observation, for all four fre- 
quencies, with a different calibrator (either PKS 1947+079 
or PKS 2029+121). For SS433, the differences in the re- 
sultant circular polarization are of the same order as the 
errors expected from system noise alone; these differences 
have been incorporated into our circular polarization error 
estimates. 

All calibrators were imaged in circular polarization as a 
consistency check, and the results for 20 May, 1999, are 
shown in Fig 1. The leakage calibration of these cali- 
brators was constrained so that the results gave zero cir- 
cular polarization for PKS 1934-638 (an absolute circu- 
lar polarization calibrator is needed for the interferome- 
ter see Sault, Hamakcr & Bregman 1996). A bias in 
the observed circular polarization of compact steep spec- 
trum sources observed at the ATCA (Rayner et al. 1999) 
suggests that PKS 1934-638 is in fact circularly-polarized 
with V/I ~ +2.5 x 10~ 4 ±0.5 x 10~ 4 at 4.8 GHz. An error 
of this order in the absolute circular polariation flux-scale 
is much smaller than the observed circular polarization 
of SS433, and does not affect the conclusions of this pa- 
per. The sign convention of V follows the IAU convention 
(Transactions of the IAU, vol 15B, (1973) 166), which con- 
forms to the IEEE definition (1969, Standard Definitions 
of terms for radio wave propogation, IEEE Trans AP-17, 
270). 

SS 433 was slightly resolved (i.e. ;> 1 arcsec) in the E- 
W direction in both total intensity and linear polarization, 



consistent with the observations of Hjellming & Johnston 
(1981). However the source is unresolved in Stokes V, con- 
sistent with all the circularly polarised flux arising in an 
unresolved point source. 

On 1999 May 10 we also observed the radio-jet X-ray 
binary GRS 1915+105; at the time the source had a flat 
spectrum between 4.80 - 8.64 GHz, at a flux density of 
~ 20 mJy. We can place 3cr limits of 1.2% on both the 
fractional linear and circular polarization of the emission 
from this source. 

In addition to our ATCA observations, we also utilise 
data from the Green Bank Interferometer (GBI) variable 
source monitoring program. These data provide a long- 
term view of the state of activity of the radio source with 
daily flux density measurements at 2.3 and 8.3 GHz (see 
e.g. Waltman et al. 1991). 

3. EVOLUTION OF THE PROPAGATING EJECTA 

Both the GBI monitoring data and our ATCA obser- 
vations are presented in Figure 2. Note that the lower 
two panels show the flux density measured in linear and 
circular polarization respectively (i.e. not the fractional 
spectra). A best fit to the observed flux density spec- 
trum of the circular polarization has a spectral index of 
ay = —0.9 ±0.1 (where V oc v av ). Between the two 
epochs it is clear that the total radio flux dropped by 
> 40% and was dominated by the decaying stages of the 
two major flare events which peaked, at 2 GHz, on ~ MJD 
51304 and ~ MJD 51312 respectively. Over the same pe- 
riod the linearly polarised flux may have increased slightly, 
and the circularly polarised flux appears to have remained 
constant. 

Our current understanding of the evolution of external 
and internal opacity effects in the outflow from SS 433 is 
summarised in Fig 3 (based upon Hjellming & Johnston 
1981, Vermeulen 1992 and Paragi et al. 1999). Also in- 
dicated are our best estimates of the physical locations of 
the components corresponding to the two flare events, at 
the times of our observations on 1999 May 10 and May 
20. We can presume that on our first epoch we only ob- 
served linear polarization from the first ejection, which 
contributed < 1/3 of the total flux density at this time. 
By our second epoch of observations both components con- 
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tributed to the linearly polarised flux. Qualitatively this 
can explain the increase in the linearly polarised flux den- 
sity between the two epochs, even though the total flux 
density decreased. Qualitatively similar behaviour is seen 
in the evolution of linear polarization in ejections from 
GRS 1915+105 (Fender et al. 1999) where the core also 
remains persistently (linearly) depolarised. 

4. THE ORIGIN OF THE CIRCULAR POLARIZATION 

Circular polarization may be produced in a synchrotron- 
emitting plasma cither directly as a result of the syn- 
chrotron process or via conversion of linear to circu- 
lar polarization (Kennett & Melrose 1998 and references 
therein). Below we briefly discuss the possible interpre- 
tations of the circularly polarised radio emission from SS 
433 in the context of these models. 

4.1. Intrinsically circularly polarised synchrotron 
emission 

An electron of Lorentz factor 7 can be considered to 
radiate synchrotron emission primarily at a frequency 
v = 4.2B±j 2 MHz, where B± is the component of the 
magnetic field perpendicular to the line of sight, mea- 
sured in Gauss. The fractional circular polarization m c (= 
Stokes |V|/J) produced intrinsically by synchrotron radi- 
ation is of order I/7 (Legg & Westfold 1968), and hence 
the observed circular polarization spectrum should follow 
the relation m c oc v^ 1 ! 2 . For an observed m c ~ 0.005 at 
8640 MHz, we can estimate a magnetic field strength of 
~ 50 mG (corresponding to 7 <~ 230). 

This estimate is in order-of-magnitude agreement with 
magnetic field estimates for major ejections from SS 433 
and other X-ray binaries (e.g. Hjellming & Han 1995, 
Fender et al. 1999). Note in addition that 'Faraday depo- 
larization', which can severely reduce the observed linear 
polarization, will have no effect on circular polarization. It 
is therefore possible to have an optically thin synchrotron 
source which displays a large ratio of circular to linear po- 
larization, as is observed in this case. However, reversals 
in the line-of-sight component of the magnetic field which 
are likely to occur within the source will require signifi- 
cantly higher magnetic fields. Similarly, a significant e + 
population within the source will also cause a reduction in 
the observed V, and require a significantly higher magnetic 
field. In addition, the ~ constant circularly polarised flux 
density at 8640 MHz during a drop by > 40% in the total 
flux density does not seem compatible with an origin for 
the circular polarization in the optically thin ejecta which 
correspond to the two major flares. 

4.2. Propagation-induced circular polarization 

Linearly polarised radiation can be converted to cir- 
cularly polarised radiation during propagation through a 
plasma with elliptical (or linear) propagation modes (Pa- 
cholczyk 1973; Kennett & Melrose 1998 and references 
therein) . In the event of the admixture of a small amount 
of relativistic plasma to a thermal plasma, propagation 
modes through the plasma will become slightly elliptical, 
and a spectrum of the form m c oc v^ 1 is predicted (Pachol- 
czyk 1973). In the event of plasma which is dominated by 
highly relativistic particles, the propagation modes may 
approach linear, and a much steeper spectrum of the form 
m c oc i/~ 3 is predicted (Kennett & Melrose 1998). 



Wardle et al. (1998) argued that the circular po- 
larization observed from 3C 279 arose because of such 
propagation- induced 'repolarization'. They concluded 
that the low-energy spectrum of the relativistic particles 
must extend to 7 << 100, and therefore the jet must be 
composed of an c + c~ plasma (if there were protons ac- 
companying each emitting electron the kinetic energy of 
the jet would be several orders of magnitude greater than 
that which is seen to be dissipated at the head of the jet). 
We note that analogous considerations may be also appli- 
cable to SS 433, where X-ray and radio hotspots are ob- 
served within the W 50 radio nebula, presumably at the 
site of the jet : ISM interaction. Indeed, the kinetic en- 
ergy in the jets of SS 433, if they are composed of a c~p+ 
plasma, is > 10 40 erg s _1 (Brinkmann et al. 1988), which 
is much greater than that which is directly observed to be 
dissipated on larger scales within W50 (this is one of the 
arguments against an c~p+ plasma which is put forward 
by Kundt 1998). 

4.3. Alternative mechanisms ? 

Given the complexity of the SS 433 system and the high 
densities and magnetic field strengths likely to be present 
near to the base of the jet, alternative origins for the cir- 
cularly polarised emission cannot be ruled out. These in- 
clude gyrosynchrotron emission from low-energy electrons 
and cyclotron maser emission (Dulk 1985). We note that 
if the circularly polarised emission is associated with a re- 
gion on the scale of one of the binary components (e.g. 

< 10 12 cm) then the brightness temperature is > 10 10 K 
at 8.64 GHz, and > 10 12 K at 1.4 GHz. 

5. DISCUSSION 

It is of great importance to determine in which of the 
various emitting regions the circularly polarised flux den- 
sity originated, in order to determine both the spectrum 
and relative strength of the emission. Our ATCA obser- 
vations rule out an extension with the optically thin jets 
on scales of > 1 arcsec, and the lack of correlated vari- 
ability in Stokes I and V argues against an association 
with the two major ejection events. The underlying ra- 
dio spectrum of SS 433 is typically around 250 mJy at 8 
GHz, with an optically thin spectral index of about —0.7, 
corresponding to a quasi-continuous flow of matter into 
the jets. If associated with this component then the spec- 
trum of the relative circular polarization could be as flat as 
the m c oc v^ 1 / 2 predicted for intrinsic synchrotron emis- 
sion. However, if associated with this component, then 
why not with the two flares, which are presumably just 
enhancements of the same flow ? Alternatively, Paragi et 
al. (1999) show that the innermost regions (< 50 mas) 
of the jets have a flat/inverted spectrum between 1-15 
GHz. The core region has a peak flux density typically of 

< 100 mJy on VLBI scales. If associated with this region 
the fractional circular polarization may be as high as 10% 
and the spectrum may steepen to the m c oc v~ x predicted 
for a mildly relativistic plasma (it seems unlikely that the 
m c oc spectrum can be recovered, unless the emission 
arises right in the binary core of the system, which has 
the most inverted radio spectrum). If associated with the 
inner regions, this implies a very large ratio of circular to 
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linear polarizations, as found for Sgr A* by Bower et al. 
(1999). 

Further precessional-phase- resolved monitoring and spa- 
tial resolution of the regions responsible for the circular 
polarization are essential to further investigate this dis- 
covery. If, as seems likely, the circular polarization is as- 
sociated with the synchrotron emitting ejecta, comparison 
with circular polarization measurements of other X-ray bi- 
naries has the potential to reveal, finally, the composition 
of the relativistic plasmas. 
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